Unconventional phonon-mediated superconductivity in MgB2 
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We have evaluated the total carrier mass enhancement factor ft for MgB2 from two independent 
experiments (specific heat and upper critical field). These experiments consistently show that ft 
— 3.1±0.1. The unusually large ft is incompatible with the measured reduced gap {2A{0) /ksTc — 
4.1) and the total isotope-effect exponent {a — 0.28±0.04) within the conventional phonon-mediated 
model. We propose an unconventional phonon-mediated mechanism, which is able to quantitatively 
explain the values of Tc, ft, a, and the reduced energy gap in a consistent way. 



The recent discovery of superconductivity near 40 K 
in MgB2 0] has led to a remarkable excitement in the 
solid-state physics community. Such a high- Tc supercon- 
ductivity in this simple intermetallic compound immedi- 
ately raises a question of whether mechanisms other than 
the conventional electron-phonon interaction are respon- 
sible for the superconductivity. A significant boron iso- 
tope effect (as = -dlnTc/dlnAlB = 0.26±0.03, where 
A-Ib is the mass of boron) and nearly zero magne- 

sium isotope effect {aMg = 0.02±0.01) 1^ suggest that 
electron-phonon coupling plays an important role in the 
pairing mechanism. To explain the 40 K superconduc- 
tivity, an electron-phonon coupling constant of about 1 
is needed. First-principle calculations give the coupling 
constant of 0.7-0.9 (Refs. I^-^), which appears to be suf- 
ficient to explain the observed high- Tc superconductivity. 
On the other hand, recent specific heat data |^ indicate 
a very strong electron-phonon coupling with a total mass 
enhancement factor of about 3.2. This would suggest a 
coupling constant of 2.2 and 2A(0)/fcsTc « 5 within the 
conventional phonon-mediated mechanism [ pO[ , in con- 
tradiction with the bulk-sensitive Raman scattering ex- 
periments which show 2A(0)//fcBTc = 4.1 (Ref. |ll|). 

Here we evaluate the carrier mass enhancement factor 
ft from the measured upper critical field. The deduced ft 
is 3.1±0.1, in remarkably good agreement with that de- 
duced from the independent specific heat data [^. The 
large mass enhancement factor is not compatible with the 
measured gap amplitude and the reduced total isotope- 
effect exponent within the conventional phonon-mediated 
model. We thus propose an unconventional phonon- 
mediated mechanism where both the long-range Frohlich- 
type and short-range electron-lattice interactions are con- 
sidered and treated distinctively. Within this scenario, 
we are able to quantitatively explain the Tc value, the 
mass enhancement factor, the total isotope-effect expo- 
nent, and the reduced energy gap. 

For a clean superconductor, as it is the case for MgB2 
the zero-temperature coherence length ^(0) is 
equal to the BCS coherence length |l^, that is. 



m 



flVp 

7rA(0)' 



(1) 



Here is the renormalized Fermi velocity, i.e., vp 
Vp/ft {vp is the bare Fermi velocity). Then 
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Using i?c2(0) = $o/27r[^(0)]2 and i?c2(0) = 15±2 
T (Ref. Il3|), we obtain ^(0) = 47±2 A. With = 
4.85x10'^ m/s (Ref. [||,0), A(0) = 6.9 meV [|ll|, and 
C(0) = 47±2 A, we find from Eq. |that ft = 3.1±0.1. Re- 
markably, the deduced mass enhancement factor ft here 
is in excellent agreement with that found from specific 
heat data 

It is clear that the large mass enhancement factor is not 
compatible with the measured energy gap (2A(0)/A:bT'c 
= 4.1) |jll[ and the reduced isotope exponent {a = 
0.28±0.04) iQ] within the conventional phonon-mediated 
model. In order to resolve the above controversy, one 
needs to consider a long-range Frohlich-type electron- 
phonon interaction that results from the Coulomb in- 
teractions of electronic charge carriers with the ions of 
ionic (polar) materials |l^]. This interaction is distinct 
from the short-range electron-lattice interaction (e.g., 
deformation-potential interaction, Holstein interaction) 
that has a carrier's energy only depending on the posi- 
tions of the atoms with which it overlaps. The Frohlich 
electron-phonon interaction is particularly strong in an 
ionic solid or a perovskite oxide where its static dielectric 
constant eo is always much larger than its high-frequency 
dielectric constant Coo , whereas the short-range electron- 
lattice interaction (including the electron-lattice interac- 
tion with acoustic phonons) is present in all the materials 
p^ . Recent Quantum Monto Carlo simulation ||l^ has 
shown that the Frohlich-type electron-phonon interaction 
is nonretarded and can always lead to a mass enhance- 
ment factor of fp = exp((7^) no matter how weak this 
interaction is. Here = A/oJo, A is a constant, and 
LOo is the characteristic frequency of optical phonons. It 
is apparent that this mass enhancement factor fp will 
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strongly depend on the isotope mass if is substan- 
tial, in contrast to the conventional retarded electron- 
phonon coupling model where the mass enhancement fac- 
tor is essentially isotope-mass independent On the 
other hand, the short-range component of the electron- 
lattice interaction is retarded and can be modeled within 
the conventional Eliashberg equations when the coupling 
constant is less than 1 and the phonon energies are much 
lower than the bare hopping integral . When only 

a short-range electron-lattice interaction is present and 
the coupling constant is far larger than 1, the interaction 
becomes nonretarded and small polarons can be formed 
[^,0. In this case, the Migdal adiabatic approxima- 
tion breaks down and one cannot use the Eliashberg 
equations to describe superconductivity |p^ , p^ . 

Since the Frohlich-type electron-phonon interaction is 
nonretarded, we can assume that the role of this inter- 
action is simply to enhance the density of states 
and reduce the direct Coulomb interaction between two 
carriers 15 1^]. In other words, we can model the re- 
tarded short-range electron-phonon coupling within the 
conventional Eliashberg equations, but the effective den- 
sity of states and the coupling constant A are enhanced 
by a factor of fp due to the presence of the Frohlich- 
type electron-phonon interaction. The effective Coulomb 
pseudopotential /i* in the Eliashberg equations will not 
change significantly when the Frohlich electron-phonon 
interaction sets in. This is because this interaction en- 
hances the density of states, which tends to increase /i*, 
while it produces an attractive potential that effectively 
reduces the direct Coulombic repulsion by a factor of 
about eo/eoo fill and thus tends to decrease ^* (see Eq. 4 
below). Therefore, the Frohlich-type electron-phonon in- 
teraction can enhance superconductivity mainly through 
increasing the effective coupling constant for the retarded 
electron-phonon interaction. For ionic materials (e.g., 
MgB2), this type of electron-phonon couphng should be 
rather strong, so that the superconductivity can be en- 
hanced substantially. 

Within this simplified approach, the effective coupling 
constant for the retarded short-range electron-phonon in- 
teraction is A = Ab/p, where \b is the bare short-range re- 
tarded electron-phonon coupling constant in the absence 
of the Frohlich-type electron-phonon interaction. The 
value of Af, for MgB2 has been calculated to be 0.7-0.9 
(Ref. [^-^). The total carrier mass enhancement factor 
is then given by 



ft = /p(l + Ab/p). 
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The factor 1 + Af,/p is the mass enhancement factor due 
to the retarded electron-phonon interaction with the en- 
hanced coupling constant A = \bfp- The Coulomb pseu- 
dopotential 11* is 
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Here ^ oc Ufp {U is the effective Coulomb interaction be- 
tween two carriers, which is renormalized by the Frohlich 
interaction ); Ep — Ep/fp {Ep is the bare Fermi energy); 
LUin is the logarithmically averaged frequency, which is 
normally lower than the Debye frequency loo by a factor 
of about 1.2. Indeed, the calculated hujin is 53.8 meV 
(Ref. which is a factor of 1.2 smaller than the mea- 

sured hujD (64.3 meV) |,§. 

Since fp — exp((7^), which depends on the isotope mass 
Ip^ , the quantities A, /i*, and ft are all isotope-mass de- 
pendent. One can easily deduce that the total exponents 
of the isotope effects on A (a a — — ^dlnA/dlnMj, 
where Mj is the mass of the jth atom in the unit cell) 
and on fi* {a\ = — '^dln fi* /din Mj) are given by 
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and 



a^* = aA[l - \ii{EF/huji„)] + fJ.*iax + -). (6) 
Furthermore, from the McMillian formula pOt], 
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we can also determine the total exponent of the isotope 
effect on {a = -Y,d\nTc/dlnMj), 
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From the above equations, we can calculate fp, A, hujin, 
a\, a^* , and a as a function of A;, using fixed parameters 

= 3.2, = 40 K, Ep = 0.57 eV (Ref. |l|]), and /i* 
= 0.1. In Fig. 1, we plot the calculated a, tiujin, and fp 
as a function of Afc. One can see that the calculated a is 
0.3 for Ab = 0.8, a value lying within the first-principle 
calculations (Afc — 0.7-0.9) The calculated isotope 

exponent a is in quantitative agreement with the mea- 
sured one (0.28±0.04) Q. The reduction in the isotope 
exponent is due to the fact that the coupling constant 
A has a negative isotope effect which partially cancels 
out the positive isotope effect on the prefactor of Eq. 7. 
As Xb increases, fp must decrease to keep ft a constant 
(see Eq. 3 and Fig. Ic). The decrease of fp reduces the 
magnitude of the negative isotope effect on the coupling 
constant, and thus increase the isotope effect on Tc, as 
seen clearly in Fig. la. 

Meanwhile, the calculated hujin is about 40 meV for Af, 
= 0.8 (see Fig. lb), which appears to be lower than the 
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one (53.8 meV) determined from the phonon density of 
states 1 19]. However, inelastic neutron scattering experi- 
ments [20 1 show that a low-energy phonon mode at about 
17 meV is strongly coupled to conduction electrons; the 
intensity of the 17 meV peak increases with decreasing 
temperature for T > Tc whereas below Tc it starts de- 
creasing. The strong coupling to the low energy mode is 
not expected from the theoretical calculations g-^], and 
may be related to some kind of structural instability . 
If this is true, the calculated hujin in Ref. [|l^ should be 
substantially overestimated. Indeed, from the measured 
2A{0)/kBTc = 4.1, one can determine the magnitude of 
hojin using a formula 
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Substituting = 40 K and 2A{0) /ksTc = 4.1 into Eq.|, 
we get hujin = 40 meV, in quantitative agreement with 
the above independent calculation. 

From Fig. Ic, one can see that fp is about 1.5 at Afc 
= 0.8, that is, the Frohlich-type electron-phonon inter- 
action enhances the density of states and the coupling 
constant by a factor of 1.5. Without this interaction, the 
material would have a transition temperature of about 22 
K. Therefore, the Frohlich-type electron-phonon interac- 
tion in ionic materials can indeed enhance superconduc- 
tivity. 

Although the conventional phonon-mediated model 
could account for the observed value 1^-^, it cannot 
consistently explain the value, the reduced isotope ex- 
ponent (~0.3), the large carrier mass enhancement factor 
(3.1), and the reduced energy gap {2A{0)/kBTc = 4.1). 
The accurate determination of 2A{0) / ksTc — 4.1 places 
a strong constraint on huJin, that is, huiin ~ 40 meV. With 
tiujin = 40 meV, and Tc = 40 K, one cannot find any pa- 
rameters that could lead to a ~ 0.28. Further, one would 
never get a carrier mass enhancement factor larger than 
3 within the conventional model. On the other hand, 
the proposed unconventional phonon-mediated mecha- 
nism naturally resolves these controversies and is able 
to quantitatively explain these experiments in a consis- 
tent way. In addition, an important prediction of this 
scenario is that both electronic specific heat 7 and Lon- 
don penetration depth Al(0) depend on the isotope mass 
of boron. One can easily show that, upon replacing 
with ^^B, both 7 and A|(0) will decrease by about 4%. 
The isotope effect on 7 should be observable if one could 
measure the specific heat of the isotope samples down to 
a low temperature (< 1 K) under a magnetic field higher 
than Hc2{0)- On the other hand, a special cation must 
be taken in determination of the intrinsic London pene- 
tration depth since the extrinsic contribution to the mea- 
sured penetration depth due to defects [ pl| would mimic 
the isotope effect if two isotope samples have different 



densities of defects. 

In summary, the carrier mass enhancement factor ft 
has been determined for MgB2 from the measured upper 
critical field and the calculated bare Fermi velocity. It 
is found that ft is 3.1±0.1, in remarkably good agree- 
ment with that deduced from the independent specific 
heat data. The unusually large ft is inconsistent with the 
measured reduced gap and the total isotope-effect expo- 
nent [a ~ 0.3) within the conventional phonon-mediated 
model. We thus propose an unconventional phonon- 
mediated mechanism where long-range Frohlich electron- 
phonon interaction and short-range retarded electron- 
phonon interaction are modeled separately. Within this 
scenario, we are able to quantitatively explain the values 
of Te, ft, a and 2A{0)/kBTc. 
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FIG. 1. The total isotope exponent a, the logarithmically 
averaged frequency Tiujin, and the mass enhancement factor fp 
(due to the Frohlich electron-phonon interaction) as a func- 
tion of the bare retarded electron-phonon coupling constant 
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